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We report an electronic magnetization measurement of a quantum point contact (QPC) based 
on nuclear magnetic resonance (NMR) spectroscopy. We find that NMR signals can be detected by 
measuring the QPC conductance under in-plane magnetic fields. This makes it possible to measure, 
from Knight shifts of the NMR spectra, the electronic magnetization of a QPC containing only a 
few electron spins. The magnetization changes smoothly with the QPC potential barrier height and 
peaks at the conductance plateau of 0.5 x 2e^//i. The observed features are well captured by a 
model calculation assuming a smooth potential barrier, supporting a no bound state origin of the 
0.7 structure. 


Quantum point contact (QPC) is a short one¬ 
dimensional (ID) channel connecting two electron reser¬ 
voirs. Its conductance is quantized to integer multiples 
of 2e^//ij_^where e is electron charge and h is Planck’s 
constant [l|, l2|. The conductance quantization is well un¬ 
derstood within a model of non-interacting electrons Q. 
However, experiments have shown an additional conduc¬ 
tance feature, a shoulder-like structure at around 0.7 x 
2e^//i termed as 0.7 structured,!^. Despite the simplic¬ 
ity of a QPC, a comprehensive understanding of the 0.7 
structure is still lacking 

Theories proposed to explain the 0.7 structure can be 
discriminated according to their predictions on the elec¬ 
tron spin arrangement, which include spontaneous smn 
polarizationd Q, antiferromagnetic Wigner crystald], 
Kondo screeningand local s pin fluctuations ac¬ 
companied by van Hove singularity [l^, [l3- Especially 
in the Kondo scenario, the existence of a localized mag¬ 
netic moment in the QPC is an inevitable ingredient. On 
one hand, early experiments observing Fano resonances 
suggested such presence of a local single spin trapped 
in a bound state regardless of magnetic fields [l^. On 
the other hand, an experiment measurin g co mpressibil- 
ity contradicts such bound state formation|18l|. Thus, the 
degree of spin polarization of a QPC is one of the central 
issues to understand the origin of the 0.7 structure. 

However, most experiments [3, IE l3“12| to date have 
focused on transmission properties, without the QPC 
spin polarization being addressed directly. Despite the 
recent progress in magnetic sensors [^, the magnetiza¬ 
tion measurement of a QPC containing only a few elec¬ 
trons is still very challenging. Recently, small magne¬ 
tizations of two-dimensional electron systems (2DESs) 
embedded in GaAs have been measured [EI-iE by com¬ 
bining tech niq ues of current-induced nuclear spin polar¬ 


ization 


2J-l28j and resistance (conductance) detection of 


nuclear magnetic resonance (NMR) signals of Ga and As 


nuclei [27H29|. Because of the hyperfine interaction be¬ 
tween electronic and nuclear spins, an electronic magne¬ 
tization produces an effective magnetic field for nuclei. 


resulting in the shift of the NMR frequency, the Knight 
shift. From the Knight shift, the electronic magnetiza¬ 
tion can be determinedfd^. 

A recent transport experiment by Ren et al.[3l[ sug¬ 
gests such influence of nuclear spins on the QPC con¬ 
ductance. They observed hysteresis in the source-drain 
voltage dependence of the differential conductance un¬ 
der magnetic fields, and attributed its origin to the dy¬ 
namical nuclear spin polarization (DNSP) induced in the 
QPC. However, NMR or other direct evidence show¬ 
ing involvement of nuclear spins has not been presented 
so far. NMR signal detection in the QPC conductance 
would constitute a novel experimental technique to probe 
spin properties of QPCs or nanowires [si [ 33 . 

In this Letter, we report an electronic magnetization 
measurement of a QPC defined in a GaAs/AlGaAs het¬ 
erostructure based on NMR spectroscopy. We find that 
the QPG differential conductance changes when the fre¬ 
quency of an applied oscillating magnetic field matches 
the NMR frequencies of ®®Ga, ^^Ga, and ^®As. The re¬ 
sistive detection of the NMR signals allows us to measure 
the electronic magnetization of the QPG from the Knight 
shifts of the NMR spectra. The Knight shift measure¬ 
ments are conducted at the QPG conductance between 0 
and 2e^//i by tuning gate and source-drain voltages. The 
magnetization changes smoothly with the QPC potential 
barrier height and peaks at the conductance plateau of 
0.5 X 2e^//i. The observed features are well captured by 
a model calculation assuming a smooth potential barrier 
without a bound state formed. Apart from the demon¬ 
stration of a new technique to measure a magnetization 
of only a few electrons, the absence of a bound state in 
the QPC is our main conclusion, directly relevant for the 
understanding of the 0.7 structure. 

QPCs studied in this work are fabricated from a 
wafer of GaAs/Alo.aGao.rAs single heterostructure with 
a 2DES at the interface. The mobility and sheet car¬ 
rier density of the 2DES at 4.2 K are 110 m^/Vs and 
2.2 X 10^® m“^, respectively. A QPC is defined electro¬ 
statically applying negative voltages {Vgi, Vg 2 ) to a pair 
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FIG. 1. (color online) (a) Linear conductance Go as a func¬ 
tion of Vgi {Vg 2 = —1.4 V) at B = 0, 1.5, 3, 4.5, and 6 T, 
applied along the x direction. Inset shows a scanning electron 
microscope image of the device, (b) Differential conductance 
G as a function of source-drain bias voltage Kd at l^i = 
—0.685 V [indicated by a dashed line in (a)] under the same 
magnetic fields as in (a). The solid (dashed) curves are mea¬ 
sured by scanning l^d in the positive (negative) direction at a 
rate of 5.6 /rV/s. Inset shows the time dependence of G at B 
= 4.5 T and Vgi = —0.685 V after an instantaneous change 
of I4d from 0 to —50 /rV. A slightly different value compared 
to the one in (b) for the same parameters, B = 4.5 T and I4d 
= —50 /rV, arises due to a remaining DNSP created at large 
|I4d| during the I4d scan in (b). 


of Au/Ti gate electrodes patterned on the surface of the 
wafer. All data presented here are measured on a QPC 
with lithographic dimensions of 300 nm length and 250 
nm width [inset of Fig.[TJa)], in a dilution refrigerator at 
the mixing chamber temperature of 20 inK. The external 
magnetic field B is applied parallel to the 2DES plane 
along the current flowing direction [x direction in the in¬ 
set of Fig. 1(a)] to avoid orbital effects and quantum Hall 
edge channels. The differential conductance G = dl/dVsd 
(where / is the current and Hsd is the source-drain bias 
voltage) is measured using a standard lock-in technique 
with a typical excitation voltage of 20 at 118 Hz. 
A single-turn coil is wound around the device to apply 
radio-frequency oscillating magnetic field Hrf. 

The QPC shows a typical conductance quantization 
behavior. Figure [TJa) shows linear conductance Gq = 
G(V[ 5 d = 0) as a function of gate voltage Vgi. In addi- 
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FIG. 2. Differential conductance G as a function of frequency 
/ of Brf at H = 4.5 T, Vgi = —0.685 V, and Fsd = —50 /rV. 
Brf is applied perpendicular to B [y direction in the inset of 
Fig. 1(a)]. / is scanned at a rate of 0.128 kHz/s. Data of 10 
subsequent measurements are averaged to improve the signal 
to noise ratio. 


tion to quantized conductance plateaus, the 0.7 struc¬ 
ture is observed at zero magnetic field, developing into 
a plateau of 0.5 x 2e^/h at high magnetic fields. A 
zero-bias conductance peak accompanying the 0.7 struc¬ 
ture is observed clearly in the G-t4d curve at i? = 0 T 
[Fig. mb)]. With increasing H, the zero-bias conductance 
peak is suppressed and turns into a dip above i? = 3 T. 

Hysteresis is observed in the G-t4d curves when I4d is 
scanned slowly (5.6 fiV/s) in the positive and negative 
directions [Fig. Iljb)]- The hysteresis is seen only at fi¬ 
nite magnetic fields. Typical time scale to develop the 
hysteresis is measured at H = 4.5 T by recording G after 
an instantaneous change of t4d from 0 to —50 /iV [Inset 
of Fig. (Ub)]- The value of G continues to change over a 
period of 200 s. This time scale is consistent with nuclear 
spin relaxation or polarization times reported in GaAs- 
based devices 27, 2^m> 34) Hi- Similarly as concluded 
in the earlier work|31j7we interpret the slow change in G 
as the first indication for the DNSP in the QPC. 

To conhrm the nuclear spin origin of the observed slow 
change in G, we perform the NMR spectroscopy exper¬ 
iment. Scanning the frequency / of i?rf, we observe 
decreases in G when / matches the NMR frequency of 
^®As (gyromagnetic ratio 7 = 45.82 rad MHz/T) [Fig. [5]. 
The obtained G-f curve represents the NMR spectrum of 
^®As, split into three dips due to the electric quadrupole 
interactionjs^. We observe signals at resonances of ®®Ga 
and ^^Ga, as well as analogous behavior in four other 
QPC devices (not shown). These observations clearly 
show that the DNSP is induced in the QPC and that its 
changes are measured by monitoring the QPC conduc¬ 
tance. 

Having established the method to probe the NMR 
spectra in transport, we now use it to determine the elec¬ 
tronic magnetization of the QPC. To this end, we per¬ 
form the following pump-probe experiment [Fig. [3][a)]. 
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FIG. 3. (color online) (a) Schematic sequence for the pump-probe experiment, (b) NMR spectra of ^®As for = 0 

fAf and various gate voltages as indicated. NMR signals corresponding to the transition between nuclear spin states 

\Iz = ±1/2) are shown. Solid curves are the fitting results. Data are offset vertically for clarity, (c) Knight shift K plotted as a 
function of gate voltage Vgi. Linear conductance Go is plotted by a dotted curve referring to the right axis. Inset shows fitting 
results of the NMR data for = —0.679 V assuming 3D (dashed) and 2D (solid) hard-wall confinement potentials, (d) 

Calculated magnetization density at the QPC center mo plotted as a function of potential barrier height Vq. The red and blue 
curves depict spin densities (no,t) and (no,4,}, respectively. Calculated conductance Go is plotted by a dotted curve referring to 
the right axis. 


First, nuclear spins are initialized by inducing DNSP un¬ 
der a relatively large bias voltage 14d = —300 fiV at Vgi 
= —0.685 V. Then, 14d is set to 0 and the QPC is 
tuned to a state of interest by setting the gate voltage 
to for a period of time (22 0, during which the 

frequency of Brf is set to / for 20 s[37|. Finally, changes 
in the DNSP are read out by recording G with a small 
ac voltage excitation (20 ^V, 118 Hz) at = —0.685 V 
and I4d = 0 fiV. The observed values of G at the begin¬ 
ning of the readout step reflect how much are the nuclear 
spins depolarized by Repeating this procedure with 
different /, we obtain an NMR spectrum for a gate volt¬ 
age as shown in Fig. [3](b). 

The bottom data of Fig. [3](b) is obtained by deplet¬ 
ing electrons from the QPC during the Rrf application. 
Therefore, this spectrum is not affected by electrons, and 
has a rather sharp dip at /o = 32.755 MHz, the frequency 
corresponding to the transition between the nuclear spin 
states |/^ = ±1/2). As is increased, the NMR 

induced dips are shifted toward negative frequencies and 
broadened. These shifts are the Knight shifts due to the 


electronic magnetizations in the QPC. 

We now evaluate the magnitude of the Knight shifts 
by taking the spatial electron distribution into account. 
Extending earlier works [Hi we adopt a model 

of electrons confined in the y and z directions with 
a transverse wave function '0(y, z). The Knight shift 
for an As nucleus at position (y, z) can be written as 
SfK{y,z) = aAsmz\i^{y,z)\'^, where oas = - 2.lx 10"^^ 
kHz m^ is the hyperfine coupling coefficiently, and 
ruz = n-f — is ID electronic magnetization density 
defined as the difference in ID spin densities. We make a 
standard assumption [y that nuclear spins are depolar¬ 
ized by the rf-magnetic field according to the detuning 
from the resonance Sf = f — (/o ± SJk) with a Gaussian 
profile exp(—(5/^/27^), where /o and 7 are the NMR fre¬ 
quency and the spectrum width without the influence of 
the Knight shift, respectively. Such depolarizations in¬ 
duce the change in the electron Zeeman energy which is 
given by an integral of local nuclear spin depolarization 
multiplied by electron distribution. Since these changes 
are small, we may expand the QPC conductance, which 
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is a function of the electron Zeeman energy, and get for 
its change 

6G{f) = A J dydz exp{-6f/2'y‘^)\ilj{y,z)f, (1) 

with A an unknown proportionality coefficient. To eval¬ 
uate Eq. o, we approximate the transverse wave func¬ 
tion by the one of a two-dimensional (2D) hard-wall 
confinement, '0(?/, z) oc cos{7Ty/ Wy) cos{'!rz /Wz) with con¬ 
finement widths Wy = (65 ± 5) nm and Wz = (18 
± 3) nm[4^. The Knight shift becomes SfK(y,z:) = 
—Kcos^{Try/wy)cos^{'!rz/wz) with a parameter K pro¬ 
portional to ruz via K = —1^(0,0)p. The exper¬ 
imental data in Fig. |3Kb) are fitted to Eq. (P) using K 
and A as fitting parameters with Jq = 32.755 MHz and 
7 = 1.36 kHz determined from the data measured at the 
depletion configuration = —0.706 V). As seen in 

the figure, the agreement of the data and the model fitted 
for each curve is excellent. 

We now consider an alternative fit, assuming 
that the QPC transport occurs through a three- 
dimensionally (3D) confined electronic state 'ip{x,y,z) oc 
cos{ttx/wx) cos{7Ty /Wy) cos{ttz/wz)- A representative re¬ 
sult, using an analog of Eq. o, is given in the inset 
of Fig. EKc) and shows a much worse compatibility with 
the data. We find that such discrepancy is not sensi¬ 
tive to the confinement details. As especially well vis¬ 
ible for large Knight shifts, the data show skewed line 
shape, with steep (gentle) slopes on the low (high) fre¬ 
quency side. This is systematically reproduced by 2D 
confinement models, unlike 3D ones (see the Supplemen¬ 
tal Material 36|). 

In Fig. EKc), K is plotted as a function of Vgi. A fi¬ 
nite K emerges near the conductance onset and increases 
steeply as the conductance is increased to 0.5 x 2e^//i. 
It keeps increasing gradually with increasing Vgi even in 
the conductance plateau region of 0.5 x 2 e^//i. As Vgi is 
increased further, K turns to decrease accompanied by 
a rise of conductance from 0.5 x 2 e^//i. As a result, a 
peak in K is formed at the high gate-voltage end of the 
conductance plateau. Using the relation between K and 
TO;;, the observed maximum value K = (11.7 ± 0.5) kHz 
corresponds to to^ = (16.5 ± 4.5) x 10® m“^. 

We now show that the observed features are well re¬ 
produced by a model calculation. We model a QPC by 
a ID tight-binding Hamiltonian, 


3,cr 




( 2 ) 

Here c] ^ creates an electron with spin a (a =t, i) at the 
j-th site of the tight-binding chain which has a hopping 
amplitude t. We assume a short-range Coulomb interac¬ 
tion represented by the on-site Coulomb energy Uj. The 
potential energy and the Zeeman energy are included in 
the on-site energy, £^- 7 / 4 , = ej ± gy.BB/2, with the Bohr 


magneton and the electron 5 -factor g. We assume 
a smooth parabolic potential barrier at the QPC center 
with a height Vq and a curvature Da,. The interaction 
term is treated by a mean-held approximation neglect¬ 
ing spin huctuations. Then the mean-held spin density 
{rij^cr) is determined by a self-consistent Green’s function 
method [ 4 ^, where the on-site energy ej^a- is shifted by 
Uj{nj^g-) with (T, the opposite spin to a. We calculate 
the magnetization density prohle mj = — Uj^i) and 

the QPC conductance Gq. The values of Uj and are 
determined from the conductance measurement data[^. 

The thick solid curve in Fig.[3Kd) depicts the calculated 
magnetization density at the QPC center toq = TOj=o 
as a function of Vq, resembling the observed Vgi depen¬ 
dence of K in Fig. EKc). According to the calculation, 
the increase in toq accompanied by the emergence of the 
conductance corresponds to the increase in the number 
of up-spin electrons in the QPC. The value of toq starts 
to decrease when down-spin electrons begin to populate 
the QPC, lifting Gq from 0.5 x 2e^ jh. The gradual in¬ 
crease in Too in the 0.5 x 2e^//i plateau region is also 
reproduced. The maximum value of the calculated mag¬ 
netization density toq = 9.3 x 10® m“^ roughly agrees 
with the value determined from the Knight shift. Spin 
polarization P = (no,t ~ + ^ 07 ) reaches 70.0 

% where toq is maximal. Distribution of rrij has a bell¬ 
shaped profile and extends over a length of about 100 nm 
around the QPC center]^. 

The gradual change in toq reflects the fact that the 
local density of states is continuous at the QPC center 
unlike in a quantum dot. We therefore attribute the ob¬ 
served gradual change in K to be consistent with a QPC 
model without any bound states. This contradicts ear¬ 
lier observations claiming that a single electron spin is 
trapped in a bound state formed in the QPCQ. We 
estimate [ 3 ^ that the observed magnitude of the magneti¬ 
zation density corresponds to the total magnetic moment 
(1.65 ± 0.45) in the QPC, exceeding the single-electron- 
spin magnetic moment which a bound state can support. 
Our measurement results of the NMR line shapes, the 
gradual change of AT, and the magnetic moment values 
are consistent with a QPC model without bound states, 
such as Refs. Bm, which predicts a smooth increase 
of the magnetization without saturation upon increasing 
the magnetic field. 

In summary, we find that the NMR signals can be 
detected by measuring the QPC conductance under in¬ 
plane magnetic fields. The resistive detection makes it 
possible to measure the electronic magnetization of the 
QPC from the Knight shifts of the NMR spectra. The 
electronic magnetization changes smoothly with the gate 
voltage and peaks at the conductance plateau of 0.5 x 
2e} jh. The gate voltage dependence of the Knight shift is 
well explained by a model calculation assuming a smooth 
potential barrier, supporting a no bound state origin of 
the 0.7 structure. 
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Supplementary Material on 

Electronic magnetization of a quantum point contact 
measured by nuclear magnetic resonance 


Electric quadrupole splitting of the NMR spectra 

The NMR spectrum of ^®As shown in Fig. 2 of the main text is split into three conductance dips. Similarly the 
spectra for ®®Ga and ^^Ga constitute of three conductance dips as respectively shown in Figs. Sl(a) and (b). These 
splittings are due to the electric quadrupole interaction, the interaction between the electric quadrupole moment 
of the nuclei and the electric field gradient at the position of the nuclei. The amplitude of the electric quadrupole 
interaction, hence the amplitude of splittings, is expected to be proportional to the electric quadrupole moment of 
nuclei. The observed amplitudes of the NMR spectrum splittings are A/ = 47 kHz, 23 kHz, and 15 kHz for ^®As, 
®®Ga, and ^^Ga, respectively. The ratio of the A/ roughly agrees with the ratio of the electric quadrupole moment 
Q in literature[l]; Qtsas = 0.29 xl0“^® ni^, QesQa = 0.19 xl0“^® ni^, and Q-noa = 0-12 xl0“^® ni^. In a cubic 
symmetric bulk GaAs crystal, the electric quadrupole interaction is usually zero. We infer that the strain induced 
by the pair of Ti/Au gate electrodes breaks the cubic symmetry and induces the electric quadrupole splittings in the 
observed NMR spectra. 


Stability of external magnetic field 

The NMR spectra shown in Fig. 3(b) of the main text are taken a± B = 4.5 T. The magnetic field is produced 
by a superconducting solenoid using a persistent mode and the field is not changed throughout the series of the 
measurements. It took about 80 hours to complete the set of the NMR spectra under various gate voltages shown in 
Fig. 3(b). We took the data for = —0.685 V twice after an interval of about 50 hours and confirmed that the 

frequency shift was less than 1.5 kHz. This indicates that the magnetic field was stable enough to resolve the Knight 
shift of the QPG. 


Hyperfine coupling constant 


Hyperfine interaction between a nuclear spin R and electron spins Sj is described by a Hamiltonian 

= E ■ u (1) 

j 

Here /jq is a magnetic permeability, pb is the Bohr magneton, and px is the magnetic moment of the nucleus which 
depends on the isotope X. is the electronic wave function amplitude at the location Ri of the nucleus. In 

a semiconductor crystal, this wave function is a product of a periodic atomic wave function u{f), and an envelope 
function namely Then we can rewrite the Hamiltonian as 


Hih) = px 


E ^I^BVxl'ipjiRz 


h = PxBk ■ Ii 


(2) 


with rjx = \u{Ri)\‘^. Values of rjx depend on nuclear species as well as the crystal structure. When electrons are 
polarized, the sum is not zero and results in an effective magnetic field (Knight field) Bk acting on the nuclear spin 
li, leading to the Knight energy shift of px\B^\, equivalent to a frequency shift (5 /k = px\BK\/h. 

For ID electrons confined in the y and 2 : directions, the Knight shift has a spatial variation reflecting the envelope 
function profile. 


SfK{y,z) = ay(n|° - nj;°)|-0(y,2)p = axmm\ip{y, , (3) 

where ax oc pxVx is a nuclear-species dependent coupling coefficient, and are spin-resolved ID elec¬ 
tron densities, and is the ID magnetization density. Assuming for simplicity a hard-wall 
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FIG. SI . Response of differential conductance G to rf-magnetic field Brf plotted as a function of frequency /. Frequency 
ranges corresponding to NMR spectra for ®®Ga (a) and ^^Ga (b) are shown, respectively. The data are obtained at the same 
condition as Fig. 2 of the main text {B = 4.5 T, Rgi = —0.685 V, and t4d = —50 /rV). 


confinement in the QPC cross-section, the transversal ground state (lowest subband) wave function is tp{y, z) = 
\j2lwy\/2lwz cosijiy/Wy) cos{ttz/wz) with confinement widths Wy and Wz- Then the maximum value of the Knight 
shift K becomes 


Similarly, in the case of 2D electron system confined in the z direction, the Knight shift becomes 

K = ax{nf^-nf^) = axm2B , (5) 

\ '^Z J \ '^Z J 

using the same coefficient ax as in the ID case. Here and n'^ are spin-resolved 2D electron densities, and 
is the 2D magnetization density. Based on the Knight shift measurements in GaAs quantum wells 
with 2D electron systems, the coefficient for X = ^^Ga is reported as = —(4.5 ± 0.2)/2 x 10^^^ kHz m^ [2] 

(This is half the original value in Ref. [2] because of the difference in the definition of a). The negative sign means 
that negative frequency shift is induced by up-spin electrons. Putting this value and values for fix and r]x = 

2.562 eh/2mp, = 1.439 eh/2mp, = 2.7x10^, and = 4.5x10^ [3, 4]) into a proportionality relation 

0^75As _ / ^75As \ / 7?75 As \ 

«7iGa \7«7iGa/ \V^iGa) ’ 

the coefficient for ^®As is determined as cttsas = —2.1 x 10“^^ kHz ni^. 


Fitting results with various confinement potentials 

In the main text, we state that a 3D confinement model can not be reconciled with the measured data. Here 
we demonstrate that this conclusion is insensitive to the wave function shape, and therefore the precise confinement 
profile. To this end, we fit the NMR spectra [Fig. 3(b) of the main text] using various wave function profiles at around 
the center of the QPG. We assume that the conductance is a function of the Zeeman energy. Because the change in 
the Zeeman energy SE^ induced by the rf-magnetic field is small, we can expand the conductance to find its change 
SG{f) being proportional to The change SE^ can be written by an integral of the local Overhauser field change 
6Bx{r; /) multiplied by electron distribution |'i/’(r^P with an unknown proportionality coefficient A. Then 6G can be 
written as 


5G{f} oc 5Ez 

= ^J df5Bf^{f;f)\tp{r)\‘^ 


(7) 

( 8 ) 


dfl[f- (/o + (5/K(fO),7]|t^{7)|^ 


( 9 ) 
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FIG. S2 . Fitting results to the NMR data, from Fig. 3(b) in the main text, for = —0.679 V assuming 6 different wave 

function profiles. The wave function profiles are listed in Table I. (a) The results for (black) and (red), (b) The results 
for ^2 (black) and (red), (c) The results for (black) and (red). The yielded fitting parameters are K = 8.9, 9.2, 
9.2, 11.3, 12.1, and 13.0 kHz and A = -4.6 x 10“®, -4.0 x 10“®, -3.4 x 10“®, -1.6 x 10“®, -1.2 x 10“®, and -1.4 x 10“® for 
#1 - #6, respectively (The wave function pre-factor C is chosen so that Max{|^|^} = 1). 


Index 

Wave function profile 

Potential 

#1 

C cos{'Ky/wy)cos{Trz/wz) 

2D hard-wall 

#2 

C exp{—y^/2ly) cos{'kz/wz) 

parabolic (y) and hard-wall (z) 

#3 

C expi-y"^/2ll)A\{z/lz + no) 

parabolic (y) and triangular (z) 

#4 

C COs{tTX/W x) COs{TTy /wy) cos{ttz/wz) 

3D hard-wall 

#5 

C exp(—a3^/2Z^) exp(—j/^/2Zy) cos{'kz/wz) 

parabolic {x, y) and hard-wall (z) 

#6 

C exp{-x^/2ll) expi-y"^/2ll)A\{z/lz + no) 

parabolic (x, y) and triangular (z) 


TABLE I. Wave function profiles used for the fitting analysis and the corresponding confinement potentials. Ai( 2 :) and no 
denote the Airy function and its largest zero point, respectively. The parameters are Wx = Wy = 60 nm, Wz = 15 nm, = ly 
= 19.8 nm, and Iz = 5.5 nm. 


where Sf^ir) = is the Knight shift at the position r. A Gaussian function = exp(—( 5 /^/ 27 ^) 

describes the depolarization of nuclear spins by rf-magnetic field at a frequency detuning Sf = f — (/o -|- 5 /k) with 
/o being the NMR frequency without the influence of the Knight shift. Rewriting /Max{\xp{r)\'^} 

and using K and A as fitting parameters, the experimental data in Fig. 3(b) in the main text are fitted to Eq. (9) 
with /o = 32.755 MHz and 7 = 1.36 kHz, the same as in the main text. 

For the fitting analysis, we assume 6 different wave function profiles '0(r) listed in table I. The fitting results to the 
y^NMR _ _Q 07 g Y (jata are shown in Fig. S2. The wave function profiles #1 - #3 correspond to a 2D confinement 
potential supposing that no bound state is formed at around the QPC center. The wave function profiles #4 - 
correspond to a 3D confinement potential supposing that a bound state is formed around the QPC center. As seen 
in the figure, the correspondence is systematically much worse in the latter cases compared to the former ones. As 
a general feature, 3D confinement models result in gentle slopes on the low-frequency side and steep slopes on the 
high-frequency side regardless of the details of the wave function profiles. This trend is opposite to the experimentally 
observed spectrum shapes, especially for those with large Knight shifts, and 2D confinement models correctly gasp 
this behavior. 

To understand this robust dependence of the conductance dip shape on the wave function confinement dimension- 
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FIG. S3 . (a) The functions Fn{p) in Eq. (10). (b) The convolution [Eq. (10)] evaluated numerically for parameters as given 
in the text for mo, mi, and m 2 = mi/20 nm, oas = —5 x 10“^^ kHz m^, 7 = 1.36 kHz, Iq =5.5, 10.5, and 13 nm for n = 1, 2, 
and 3, respectively. 


ality, we rewrite the spatial integral in Eq. (9) into an integral over the electron density p = 

SG{f) oc dpFn{p)exp ( ——— ^72 ) ’ 

where we introduce the magnetization densities for two (n = 2) and three (n = 3) dimensional confinement as 
mi = 16.5irm^^ (our experiment) and mo = 1 (a polarized local moment), respectively. The distribution of the 
density p between its minimal p = 0 and maximal value p = pM is given by the weight function Fn{p). For 
illustration, we take a spherically symmetric Gaussian density profile, |\I'(r)p = pm exp(—r^/Z§), for which Fn{p) 
can be calculated analytically as Fn{p) cc [ln(pM/p)]^" This function is plotted on Fig. S3 (a), from where the 

dimensionality effect is clearly visible. Indeed, for n = 2 (red curves) the weight function is constant, so its convolution 
with the Gaussian frequency distribution will result in a symmetric profile, as shown on Fig. S3 (b). Gompared to 
that, a three dimensionally confined wave function has more weight at small p, which in other words means that 
it has a more weight in the tail. This results in an asymmetric shape of the convolution, skewed to the left (blue 
curves). Because ctas is negative, the shape is mirror reflected when converting the horizontal axis from the Knight 
shift in Fig. S3 (a) to the NMR frequency in Fig. S3 (b) and Fig. S2. This allows us to exclude the possibility of a 3D 
confinement model, as our data systematically show either symmetric, or skewed towards lower frequencies shapes, 
corresponding to 2D, and ID confinement models, respectively. 

We also note that the black curves - ^3) in Fig. S2 look very similar, suggesting that the fitting does not 
depend strongly on the details of the wave function profiles, and yield robust values of K, being 8.9 kHz, 9.2 kHz, 
and 9.3 kHz for #1, and #3, respectively. There is a recognizable trend, namely the values of K tend to increase 
as the tails of the wave function profile extends. Therefore the model of the 2D hard-wall potential we adopt in the 
main text gives a bound from below for K, hence for which already exceeds the magnetic moment for a single 
electron spin. This further corroborates that our conclusions are not artifacts of a specific confinement model. 


Model calculation 


In our tight-binding model calculation, we assume a smooth parabolic potential barrier at the center of the QPG 
with a height Vq and a curvature Intermediate regions {jo < \j\ < N) are inserted between the central part of 
the QPG and the leads, where ej = 0, to connect the potential smoothly [Fig. S4(a)][5]; 


|j|<Jo 

a{j - Nf + b{j - N)\ jo < \j\ < N 
0, N< \j\ 


( 11 ) 
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FIG. S4 . Tight-binding model parameters for N = 25, jo = 13, Qx/t =0.078, Vo/t = 0.586, and U/t = 0.93. (a) Potential 
energy as a function of site index j. (b) Site dependent Coulomb interaction Uj as a function of j. 
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FIG. S5 . (color online) (a) Transconductance (dG/dV^i) plotted in color scale as a function of t4d and Pgi for B = 0.0 T. 
The dotted lines trace the alignment of subband edges with the source and drain electrochemical potentials. A lever-arm factor 
converting from gate voltage to energy is derived from their slopes and intersections, (b) Transconductance for t4d = 0 /LtV and 
f? = 4.5 T plotted as a function of Vgi. The strength of Coulomb interaction U is set by requiring that the calculated result 
reproduces the measured energy difference between the spin-split transconductance peaks. 


where a and b are determined so that cj and dej/dj are continuous at j = jo'. 

a = (^0 - + %JoOo “ ^)) . (12) 

^ = 407^ (-fio - 2a(jo - iV)) . (13) 

The on-site Coulomb energy Uj is non-zero only in the constriction part and changes smoothly to Uj = 0 in the lead 
[Fig. S4(b)] [5, 6]. 


f U exp ( 

' iJ/Nf \ 

, i-u/Ny^J 


1 0, 


IjI > N. 


(14) 


The parameters for the model calculation are determined as follows. The hopping amplitude t is related to the 
Fermi energy ep by ep = 2t(l — cos27ra/Ap), where Ap is the Fermi wave length and a is the lattice constant of the 
tight-binding chain. The numerical calculations are done for N = 25 and jo = 13. The lattice constant a is set to 
a = l/8Ap [7], which leads to the QPC length of 340 nm for N = 25. A lever-arm factor rj = 34 meV/V converting 
gate voltage Vgi to energy is extracted from the transconductance (dG/dVgi) data in Fig. S5(a). Then, assuming 
a transmission probability t„{e) = 1/(1 -|- the curvature of the potential barrier is evaluated 

by fitting the linear conductance transition between 2e^ jh and jh. The parameter U is determined so that the 
calculated result reproduces the energy difference between the transconductance peaks corresponding to the spin-split 
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FIG. S6 . Calculated magnetization density nij as a function of site index j for (ep — V(,)/Q,x = 0.82 where mj=o is set to the 
maximal value observed in Fig. 3(d) of the main text. 


conductance transitions [Fig. S5(b)]. For the QPC measured, we find t = 12.8 meV, fla, = 1.0 meV and U/t = 0.93. 
The value of g-factor j^j = 0.44 for bulk GaAs is used for the Zeeman energy. 

Figure S6 depicts the spatial distribution of the calculated magnetization density mj for (ep — Vo)/Qx = 0.82 where 
mj=o is set to be maximal seen in Fig. 3(d) of the main text. The distribution is a bell-shaped curve at the QPC 
center with a length at half maximum of 16 sites which corresponds to the distribution length L of about 100 nm. 
We use the calculated L to convert the observed magnetization density to the total magnetic moment M in the 
QPC, by M = m^L. In this way we get the result stated in the main text, that the observed magnetization density 
of rtiz = (16.5 ± 4.5) ni“^ corresponds to a total magnetic moment of M = 1.65 ± 0.45. 
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